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ABSTRACT 

The two-year COBE-DM.H 53 and 90 GHz sky maps are used to determine the nor- 
malisation of inflationary, flat, dark matter dominated universe models. The appropriately 
normalized cold and mixed dark matter models, computed for a range of values of O5 and 
h, and several values of the hot to cold dark matter abundance ratio, are then compared 
to various measures of structure in the universe. 
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1. Introduction 

The COBE-DM.H discovery of cosmic microwave background (CMB) anisotropy 
(Smoot et al. 1992, Bennett et al. 1992, Wright et al. 1992) has affected cosmology in 
both ontological and practical ways, but its predominant quantitative influence has been 
to provide the means for the accurate normalisation of theories of large scale structure 
formation. 

The development of inflationary ideas during the 1980s induced a decade-long adher- 
ence to the cosmological paradigm which posits that the universe is spatially flat. Such 
a picture requires that the present energy density of the universe is dominated by non- 
baryonic dark matter or alternatively by a non-zero vacuum energy contribution (a cos- 
mological constant term, A). The minimal version of the model, which invokes cold dark 
matter (CDM) as the major constituent of the universe, is presently in direct confrontation 
with astronomical observations. An extension of the model, which in addition to CDM 
postulates an admixture of hot dark matter (HDM), enjoys considerable popularity in con- 
temporary cosmological research. Whilst vigorous discussion ensues in the literature as to 
the plausibility of the mixed dark matter (MDM) model as a viable cosmology, (see e.g. 
Schaefer, Schafi & Stecker 1989, Klypin et al. 1993, Pogosyan & Starobinsky 1993, Ma & 
Bertschinger 1994, Primack et al. 1994), we herein consider both CDM and MDM versions 
of the inflationary scenario (with A = 0). 

In this Letter we use the linear angular power spectrum estimation technique (Gorski 
1994) to normalise the cold and mixed dark matter models to the two year COBE-DM.H 
anisotropy measurements (Bennett et al . 1994), and subsequently discuss some predictions 
for large-scale structure measures resulting from this normalisation.! 

2. CMB Anisotropy Normalisation Procedure 
2.1 Data Selection and Power Spectrum Inference Method 

The two year COBE-DM.H 53 and 90 GHz sky maps are used identically as in Gorski 
et al. (1994). In addition to the galactic sky maps we also utilise the ecliptic coordinate 
frame data sets as a check on the extent to which the coordinate dependent noise binning 
can affect the inferred normalisation. 

We implement the power spectrum estimation method of Gorski (1994). Coordinate 
system specific orthogonal basis functions for the Fourier decomposition of the sky maps are 

f Similar considerations were recently presented in Bunn et al. (1994). 
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constructed so as to include exactly both pixelisation effects and a |6| < 20° galactic plane 
excision (leaving 4016 and 4038 pixels in the galactic and ecliptic sky maps, respectively). 
A likelihood analysis is performed as detailed in Gorski et al. (1994). 

2.2 Theoretical Spectra of Anisotropy 

The cosmological models are specified to within an arbitrary amplitude as follows: 
1) the global geometry is flat with A = and fio = 1, with the bulk mass density provided 
by either cold or mixed dark matter (CDM or MDM); for the MDM models the hot 
dark matter is introduced in the form of either one or two (equal mass) families of massive 
neutrinos, with the contributed fraction of critical density taken as fi„ = 0.15, 0.2, 0.25 and 
0.3 for one massive flavour, and fi„ = 0.2, 0.3, otherwise; 2) for the CDM model the values 
of the Hubble constant, Ho = 100 h km s _1 Mpc -1 , and baryon abundance are sampled at 
h = 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8, and fi 6 = 0.01, 0.03, 0.05, 0.07, 0.1 (for all values of h), 
respectively; in order to trace the Big-Bang nucleosynthesis (BBN) relation, O5 = 0.013 h~ 2 
(Reeves 1994), we also use 6 = 0.02 for h = 0.8, and 6 = 0.14 for h = 0.3; for the MDM 
model h = 0.5 and O5 = 0.05; 3) random-phase, Gaussian, scalar primordial curvature 
perturbations (no gravity waves) are assumed with the inflationary Harrison-Zel'dovich 
spectrum corresponding to an adiabatic density perturbation distribution, P{k) oc k. 

The CMB anisotropy multipole coefficients and the matter perturbation transfer func- 
tions for both models were evaluated using the Boltzman equation code of Stompor (1994) 
by solving the propagation equations up until the redshift z = 0. Over the low-£ range of 
CMB multipoles probed by COBE-DM.H the theoretical spectra are indistinguishable for 
CDM and MDM models with equivalent h and O5. Thus, the power spectrum amplitude 
derived from the data applies equally to both CDM and MDM models. 

Figure 1 shows a selection of the CMB power spectra normalised to the COBE- 
DMR two year CMB anisotropy. We parameterize all such spectra by the exact value 
of the quadrupole, denoted Qrms-PS , hi a straightforward generalization of the Qrms-PS 
introduced in Smoot et al. (1992) for pure power law model spectra. This exactly evaluated 
present-day quadrupole is of smaller amplitude (by ~ 10%, 5%, or 2% for h ^ 0.4, h ~ 0.5, 
and h ~ 0.8, respectively) than its pure power law, Sachs & Wolfe (1967) counterpart, 
mainly due to the high-redshift, integrated Sachs- Wolfe effect. 
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2.3 Results of Qrms-PS Fitting 

A typical likelihood fit of a flat dark matter model to the two year COBE-DM.H data 
yields a ~ 13a significant determination of Qrms-PS ~ 20 /iK. Systematic shifts in the 
central value of the fit are observed due to: 1) differences in the noise pixelisation in the 
galactic and ecliptic coordinate frames, which result in a ~ 0.8 /iK difference between the 
inferred normalisation amplitudes, with higher values obtained from the ecliptic maps; 2) 
exclusion of the quadrupole, which produces a ~ +0.4 /iK variation in the fitted amplitude; 
3) the uncertainty in the values of h and O5 (reflected by differences in spectral shape over 
the ^-range accessible to DMR) causes an additional small spread, ±0.25 /iK, of the fitted 
amplitudes. 

A convenient summation of the proposed overall normalisation for the flat, dark matter 
models is then Qrms-PS = (20±1.52)(±0.4 ±0.2 ±0.25)/iK. The error ranges represent the 
statistical error and uncertainties associated with effects 1 through 3 above, respectively. It 
will be noted that the statistical error on the inferred normalisation is considerably larger 
than the other uncertainties. 

Gorski et al. (1994) showed that for the power law models specified by Qrms—PS 
and the spectral index n a convenient, n-independent normalisation was in terms of the 
multipole amplitude 09 ~ 8 /iK. This ^-order is related to the point at which the theoretical 
signal to noise ratio is of order unity. For the dark matter models discussed here, we find 
that the appropriate pivot point is at an ~ 7.15 ± 0.55/iK. 

3. Measures of Large-Scale Structure 
Having determined the normalisations for our grid of CDM and MDM models, we 
now proceed to discuss the values of several large-scale structure statistics computed from 
the matter perturbation spectra according to the usual prescriptions given as footnotes to 
Table 1. Representative values for the h = 0.5, O5 = 0.05 models are found therein. 

3.1 Mass Fluctuations: a% and J3 

COBE-DM.H normalised values of (as)mass are shown in Fig. 2a. A related, observ- 
able quantity is the fluctuation in the number of galaxies within a sphere of fixed radius. 
Recent estimates thereof are close to the standard value from Davis & Peebles (1983), 
(as)gai — 1 — a representative selection of the galaxy surveys is: (as) ga i — 0.8 (Fisher 
et al. 1993), 0.87 ± 0.07 (Feldman, Kaiser & Peacock, 1994), and 0.83too 7 (Baugh & 
Efstathiou, 1993). 
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In order to compare the rms mass fluctuations with the rms galaxy fluctuations, we 
adopt the simple picture of biased galaxy formation, wherein ((Js)gai = b((T%) mass , and b is 
the linear bias factor. Some constraints on the value of b can be imposed by recent galaxy 
surveys, but are more uncertain than the ((Js)gai determinations. The estimated values of 
b range from 0.9 to 2, with strong bimodal behaviour about b ~ 1 and 2 (see Table 1 in 
Dekel 1994, and Table 1 of Cole, Fisher & Weinberg 1994, with fi = 1). 

Estimates of ((Js)mass inferred from galaxy cluster catalogues favour a higher value for 
b. Henry & Arnaud (1991) used the abundance of clusters as a function of X-ray temper- 
ature to derive ((Js)mass = 0.59 ± 0.02 for a scale-free, flat universe. White, Efstathiou & 
Frenk (1993) have used the masses and abundances of rich clusters of galaxies to determine 
(c r 8)mass — 0.52 - 0.62 for a critical density universe (relatively independent of h). A low 
value of ((Js)mass ~ 0.5 is also required by the observed low pair-wise velocities (Davis et 
al. 1985, but see Zurek et al. 1994). 

Although the value of the rms mass fluctuations is not yet well known, it seems 
unlikely that the biasing parameter is less than unity. This would require either pure 
CDM with a \ow-h or high-fij, in disagreement with observations (0.4 ^ h ^ 0.8) or the 
BBN constraints, respectively, or alternatively an appreciable admixture of a hot dark 
component. 

We have also computed the ^-integral on a scale R = 20/i _1 Mpc (Fig 2b). This, in 
principle, is less contaminated by non-linear galaxy power spectrum evolution than a%. The 
predictions of J3 are compared to the value Js(20h~ 1 M.pc) ~ 700 h~ 3 Mpc 3 from Davis & 
Peebles (1983). 

3.2 Galaxy and Matter Distribution Power Spectra 

In principle, the most direct and informative comparison of theory and observations 
can be conducted using the predicted and empirical power spectra. Several recently esti- 
mated galaxy power spectra are shown in Fig. 3, together with various theoretical linear 
power spectra (assuming b =1). However, such a comparison is still subject to difficulties. 
There are notable discrepancies in both shape and amplitude between galaxy power spec- 
tra derived by varying techniques from different samples. This may be due to the different 
power spectrum estimators employed, or intrinsic variations in the samples. Furthermore, 
on small scales, the effects of non-linear evolution should be apparent, and possible distor- 
tions can be introduced by redshift space mapping. 
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With these caveats, we only proceed to use Fig. 3 illustratively. It is apparent that 
none of the empirical power spectra exhibit as much small-scale power as the linear C OBE- 
DMR normalised CDM spectrum f . This could be considered as supportive of the MDM 
models. However, the high amplitude of the matter power spectra implied by the COBE- 
DMR measurements leaves very little room for the biasing parameter to exceed unity, in 
conflict with the 6-estimates from cluster properties. Baugh & Efstathiou (1993) have 
suggested that, with b ~ 2, the matter power spectrum required by the observed galaxy 
distribution should look rather like that for the = 1 } T = 0.2 model |, depicted in Fig. 
3 as a thin solid line. Alas, none of the CO-B-E-DMR-normalised CDM or MDM model 
power spectra resembles such an ad hoc spectrum. 

3.3 Large-Scale Flows 

The local streaming motions of galaxies provide an interesting constraint for cosmolog- 
ical models. In particular, galaxy peculiar velocities directly measure mass fluctuations, 
independently of a linear bias parameter. Dekel (1994) gives estimates of the average 
peculiar velocities within spheres of radius 1000 to 6000 km s _1 . These seem to be in 
agreement with the HI data from Giovanelli & Haynes (Dekel, private communication). 
An important recent development in the field was the determination by Lauer & Postman 
(1994) of significant bulk flow in a deep volume limited sample of ~ 100 galaxy clusters. If 
confirmed, this observation would likely invalidate all of the presently considered models 
of structure formation, which predict too rapid a decrease in the bulk flow amplitude with 
scale. The COBE-DM.H normalised rms bulk flows are shown in Fig. 2c. The model 

f This conclusion is based on the real space power spectrum estimates from Baugh 
& Efstathiou (1993) and Peacock & Dodds (1994), and the redshift space power spectra 
derived from the IRAS 1.2 Jy survey. However, the redshift space spectra from da Costa 
et al. (1994) indicate yet more power on small scales, a discrepancy which can only become 
more pronounced when one recalls that the real-to-redshift space mapping suppresses power 
on scales k £ 0.1/iMpc -1 (G ramann et al. 1993). Similarly, the predicted linear real-to- 
redshift space power enhancement on scales k ^ 0.03/iMpc -1 ) is too small to account for 

the amplitude discrepancies in these data sets for b ^ 0.4. 

| Here, we follow the parametrization of the matter power spectra introduced by Efs- 
tathiou, Bond, & White 1991. 
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predictions are in good agreement with the POTENT data, while the disagreement with 
the Lauer & Postman result appears rather significant. 

4. Discussion 

The improved quality of the two year C OBE-DM.H data combined with reliable power 
spectrum estimation techniques allows the accurate normalisation of cosmological theories. 
Previously, the variance of COBE-DM.H temperature fluctuations on a 10° angular scale 
was utilised for the normalisation of the power spectrum. Subsequent work (Wright et al. 
1994, Banday et al. 1994) has demonstrated that this technique can be unreliable without 
considerable attention. More appropriate methods take full advantage of the measured 
CMB anisotropy power distribution on all angular scales accessible to the COBE-DM.H 
instrument, as implemented in this Letter. 

The cold dark matter theory with a standard choice of cosmological parameters re- 
quires a very high normalisation in order to fit the CMB anisotropy distribution. Analysis 
of the first year of COBE-DM.H data had already suggested that ((Js)mass ~ 1 (Wright 
et al. 1992, Efstathiou, Bond, & White 1992), and this value increases to ~ 1.4 with two 
years of data and an improved analysis technique. Although this normalisation allows the 
theory to predict large-scale velocities of comfortably high amplitude, it also results in a 
significant overproduction of density perturbations on scales of ^ 20/i _1 Mpc. CDM has 
often been criticised for its poor match to both galaxy and cluster distributions. Mixed 
dark matter models manage to circumvent, to a certain degree, these same problems by 
construction — massive neutrinos partially damp the perturbations at those length-scales 
where CDM looks problematic. Among the MDM models those with two species of massive 
neutrino seem to meet the observational constraints more comfortably (see also Primack et 
al. 1994). The larger free-streaming radius allows for the suppression of the perturbation 
amplitude on larger scales than in models with one massive flavour. This is reflected in 
the decrease of the predicted ((Js)mass values. Nevertheless, the proponents of MDM will 
have to address the viability of the model viz. the simultaneous requirements that there 
should be no bias between the galaxy and mass distribution (as suggested by this analysis) 
and the galaxy pair-wise velocities should be small. 

Whilst it would be safer to await the final 4-year COBE results before offering defini- 
tive statements as to the viability of theoretical models, one should note that the CDM 
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normalisation derived from the two year COBE-DM.H data does appear to be irreconcil- 
ably high, while the MDM model has little room left for adjustment. 

We acknowledge the efforts of those contributing to the COBE-DM.H. This work was 
supported in part by the Office of Space Sciences of NASA Headquarters. We are grateful 
to L. da Costa, K. Fisher, M. Vogeley for providing their power spectra, and A. Dekel for 
providing the bulk flow points. 
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FIGURE CAPTIONS 

Fig. 1. — CDM/MDM CMB anisotropy power spectra normalised to COBE-DM.H (where 
Af = £(2£ + l)a|/47r). The four sets of curves correspond to the power spectrum ampli- 
tude determination as follows (from top to bottom): 1) ecliptic coordinates, quadrupole 
excluded, 2) ecliptic coordinates, quadrupole included, 3) galactic coordinates, quadrupole 
excluded, and 4) galactic coordinates, quadrupole included. Within each set the three 
curves denote: h = 0.5 and O5 = 0.05 - heavy type, h = 0.3 and O5 = 0.14 - medium- 
heavy type, h = 0.8 and O5 = 0.02 - light type, which are consistent with the constraints 
from BBN. The crossing point within each set is at £ 11. Several pure Sachs- Wolfe power 
law spectra are shown (shifted down to 0.9 at £ = 2). Although a power law approximation, 
P(k) oc k n , used to generate the multipole coefficients af solely through the Sachs & Wolfe 
(1967) effect will be poor for such spectra, a value of n ^ 1.1 would be most appropriate 
over the range £ ^ 15. This is a little steeper than the underlying, inflationary n = 1 
spectrum (Bond 1993). The proposed overall two year COBE-DM.H normalisation (trans- 
lated from Qrms-PS ~ 20 /iK) is represented by a filled circle; the error bar represents 
a typical statistical uncertainty of an individual likelihood fit. Clearly, this encompasses 
the fit uncertainties due to choice of coordinate system, quadrupole inclusion/exclusion, 
and/or cosmological parameter values. 

Fig. 2. — a ) ((Js)mass values predicted from the C OBE-DM.H normalised, flat, dark matter 
dominated models. The thin solid lines correspond to the CDM model with h decreasing 
from 0.8 to 0.3 in steps of 0.1 from top to bottom. The thick solid line shows the models 
which obey the Big Bang nucleosynthesis constraint, O5 = 0.013 h~ 2 . The individual points 
represent several mixed dark matter models: triangle - fi„ = 0.15, square - fi„ = 0.20, 
hexagon - fi„ = 0.25, diamond - fi„ = 0.30, where the filled symbols have N v = 1, open 
symbols N v = 2. All of the points correspond to h = 0.5 and O5 = 0.05 but are spread 
out on the plot for clarity. 

b ) The values for the J3 integral over the density perturbation correlation function within 
20 /i _1 Mpc (units [/i _1 Mpc] 3 ). Same coding as a. 

c) rms amplitudes of the large scale flows ([kms -1 ]). The heavy lines correspond to the 
CDM models with h = 0.8, fi 6 = 0.02 — top, h = 0.5, fi 6 = 0.05 — middle, h = 0.3, 
O5 = 0.14 — bottom. Thin lines separating from the h = 0.5 curve correspond to the 
MDM models with N v = 1 (higher line), and N v = 2, practically independent of fi„. 
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POTENT data (courtesy A. Dekel) are shown by circles, and the square shows the Lauer 
& Postman datum. 

Fig. 3. — COBE-DM.H normalised inhomogeneity power spectra and, shown as the lower- 
most curves, miscellaneous spectral windows required for the computation of the statistics 
considered in this paper (as annotated). The theoretical mass distribution power spectra 
for the fio = 1, h = 0.5, O5 = 0.05 models are shown in the middle section of the plot: 
the top heavy line shows the CDM power spectrum; the lower heavy lines correspond to 
two N v = 1 MDM models — fi„ = 0.2 (upper) and fi„ = 0.3 (lower); the medium heavy 
lines correspond to the equivalent N v = 2 MDM models. Note that the MDM transfer 
function for N v = 2 drops faster near k ~ 0.04 than for N v = 1. The thin line shows the 
r = 0.2 transfer function for CDM. The vertical error bar on the P{k) ~ k part of the 
spectrum above the an window illustrates the ~ 1 a allowed variation in amplitude for 
this fixed slope. Conversely, the superposed 'bow' shows the allowed 1 cr-variation (~ 0.3) 
in tilt for spectra with the fixed an = 7.15 /iK amplitude. It should be noted that this is 
the proper representation of two dimensional uncertainty in the spectrum determination 
from the two year COBE-DM.H data for flat dark matter models. One is not at liberty 
to arbitrarily vary both amplitude and spectral slope simultaneously. This is equivalent 
to the observed degeneracy seen in the two-dimensional (Qrms-PS , n) fits to pure power 
law models (Gorski et al. 1994 and references therein). The upper section of the plot 
reproduces the theoretical spectra for CDM, MDM (N v = 1, fi„ = 0.2), and T = 0.2-CDM 

- all shifted upward for clarity. Overplotted are several observational estimates for the 
galaxy distribution power spectrum: squares - Baugh & Efstathiou (1993), filled circles 

- Peacock & Dodds (1994), pentagons - Fisher et al. (1993), triangles - da Costa et al. 
(1994). 
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Table 1: Inferred cosmological statistics for models with Slo = &cdm + &v + &b = 
1, Cl h = 0.05, h = 0.5, and a COBE-DMR normalisation of Q rms -ps = 20.04 fiK. N v 
is the number of massive neutrino species, and m v the neutrino mass in eV. The errors, 
including both statistical (lc) and systematic deviations, are of the order of 11%. 
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